Transport through a quantum dot with two parallel Luttinger liquid leads 
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We study transport through a quantum dot side-coupled to two parallel Luttinger liquid leads 
in the presence of a Coulombic dot-lead interaction. This geometry enables an exact treatment of 
the inter-lead Coulomb interactions. We find that for dots symmetrically disposed between the two 
leads the correlation of charge fluctuations between the two leads can lead to an enhancement of 
the current at the Coulomb-blockade edge and even to a negative differential conductance. Moving 
the dot off center or separating the wires further converts the enhancement to a suppression. 

PACS numbers: 71.10.Pm, 73.21.-2b, 73.63.Kv, 73.63.Nm 



I. INTRODUCTION 



Quantum dots and single-molecule devices reveal many 
quantum phenomena of fundamental and perhaps prac- 
tical interest.—"— However, the role of a Coulombic inter- 
action between the charges on the dot and on the leads 
has been little studied, even though this Coulomb in- 
teraction is clearly relevant. Dot-lead interactions were 
studied by Boulat and Salcur— , and Boulat, Salcur and 
Schmitteckert^ in the context of an interacting resonant- 
level model. Goldstein and collaborators presented gen- 
eral duality-based results that they noted were applicable 
also to models with dot-lead interactions In previous 
papers we considered different aspects of quantum dots 
Coulomb-coupled to leads, including the relaxational dy- 
namics of a quantum dot side-coupled to a single Lut- 
tinger liquid lead in equilibrium (in Rcf. 0) and the case 
of a quantum dot coupled to two one-dimensional leads 
subject to a finite bias voltage (in Ref. Q). However, in 
the theoretical work to date inter-\ead couplings have 
not been considered. While inter-lead interactions may 
be negligible if the screening of the Coulomb potential 
is strong, in general one expects molecular junctions to 
involve quantum dots which are very close to the leads, 
so that the inter-lead interaction is not negligible. In this 
paper we show that this interaction can be of crucial im- 
portance, because it suppresses fluctuations which act to 
decohere the different states. Thus, electron-electron in- 
teractions can enhance transport, leading (in appropriate 
circumstances which we define below) to conductances 
which peak at the Coulomb-blockade edge and thus ex- 
hibit a negative differential conductance. 

A lead may be modelled as a quantum wire described 
by Luttinger liquid physicsi^"— However, for general ge- 
ometries, inter-lead interactions do not fit easily into the 
Luttinger liquid formalism, because the boundary breaks 
translational invariance and the interactions give rise to 
a complicated boundary condition on the charge modes. 
We are not aware of a convincing theoretical treatment 
of this boundary condition. In this paper we therefore 
specialize to the geometry shown in Fig. [1] for which the 
inter-lead interaction may be treated by standard Lut- 
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FIG. 1: Schematic of a one-dimensional lead-dot-lead system 
showing two parallel leads. 



finger liquid methods. This geometry reveals the essen- 
tial physics of the inter-lead coupling. 

The rest of this paper is organized as follows. Section 
im defines the model we study and the methods we use. 
Section IIIII gives the basic theoretical results. Sec. IIVI 
presents the experimental consequences, and Sec. |V] is 
a conclusion. 



II. MODEL AND METHODS 
A. Model 

The quantum dot problem sketched in Fig. [T] is de- 
scribed by a Hamiltonian of the general form 



H = H\ 



lead 



dot 



Coul 



(1) 



We label the two leads hy a ~ L , the lead orbitals by 
a and write the lead Hamiltonian as 

a—L,R aka 

^Jdxdx' ^ Vga' {X - x')paa{x)paUi' {x'). 

aa' —L,R aa' 

(2) 



The operator (^^^^^^ creates an electron with momentum 
A:, energy e^, and spin a in state a of lead a. The inter- 
action V depends on the difference in position x~ x' and 
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is the sum of two terms, an intra-lead interaction 

^2 



(x X ) 



e\x — x'\ 



(3) 



(already considered in Ref. |8|), and an inter- lead interac- 
tion denoted by 



Vj^Q,(x X ) 



i^{x - x'Y + (P 



(4) 



with a = L \i a = R and vice versa. paa{x) is the corre- 
sponding operator giving the charge density at position 
x, e is a background dielectric constant, and d denotes the 
spacing between the leads. The convenient feature of the 
geometry we consider is that both intra- and inter-lead 
interactions are functions only of the relative distance 
along the lead, permitting use of usual bosonisation tech- 
niques. 

The quantum dot Hamiltonian -ffdot may be written as 



Hdot = £d nd + ^nd{nd - 1), 



(5) 



where U is the dot charging energy, Ud = d^d^ is the 
total dot density, and dj. creates an electron with energy 
Ed and spin a on the dot. 

The Coulomb interaction between electrons on the 
quantum dot and electrons in the leads may be written 
as 



Hcoul = nd X! X! / dxVVa{x)paa{x) 
a=L,R a •' 



with 



Wo^ix) 



Va;2 + dl 
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(7) 



and da the distance between the quantum dot and lead 
a. The dot- lead hybridization is given by 



dl 



a—L.R aka 
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B. Bosonisation 



We bosonisc the system as in our previous papers-^ 
the new feature is the inter-lead coupling, which mixes 
the boson modes of the two leads. We consider a system 



with linear dimension L and periodic boundary condi- 
tions, and combine spin and orbital quantum number 
into a superindex j3 = 1, . . . , Af. The physics is repre- 
sented in terms of right (A = and left (A = — ) moving 
particle-hole pairsji^ which can be recombined into boson 
operators (j^apil), 110,^(5) that obey the volume commu- 
tation relation [(j)ap{q),Tla'l3'i~q')] = i{L/TT)daa'Spp'Sgq'. 

The total particle density in lead a and state /3 is given 

by Papiq) iq(j)af}{q)- 

The lead electron creation operator VqA^ may be 
rewritten in terms of the bosons a.^^^ 



'[A0„^(9)-in„3(q)] 



(9) 

Here the small positive infinitesimal factor 77 arises from 
the correct normal ordering of the operators. The op- 
erator UaX denotes the Klein factor which carries the 
Fermi statistics. Left and right lead electrons are as- 
sumed to be at different chemical potentials p.^, which 
gives rise to a different time evolution of the Klein fac- 
tors, Uax{t) = e'^°*C/„A(0). 

The symmetry between the two leads means that their 
low-energy physics can be described in terms of a set 
of new boson operators ipab^^ab with a parity quantum 
number p = e, o as well as a mode label b. The new 
operators are related by a linear transformation to the 
0Q/3,na/3 defined in terms of the fermions. In terms of 
the new operators the lead Hamiltonian becomes 



i^lcad = E E E^^P''^'?) 



p=c,o t=l,...,A/ q 



Kpb{q)Ilpb{-q)Ilpb{q) 



Kpb{q) 



<l>pb{-q)(f>pb{q) 



(10) 



with even (e) and odd (o) boson modes defined by 



TT f \ n f \ IlLbiq) ±TlRbiq) 
Hoh((7),Hoh(g) = -j= , (11) 



and 



, . X , f . 4>Lbiq) ± 0flfc(g) 
(pcb[q),(Pob[q) = , (12) 

and Luttinger parameters Vcb, Vd, and Kcb, Kob deter- 
mined by the bare velocities and interactions of the lead 
eigenstates. In terms of these boson modes, the electron 
operator becomes 
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with 'i/'^A/j s-n exponential of a combination of the 4>a,b, 
IIq b with b > 2 

Rewriting the Coulomb potential as 



(15) 



with A a short-distance cutoff and 
1 



Y E e''" (9) , 14a (g) = log ( 1 + 



we find for the Luttinger parameters in Eq. ([TU)) that 

1 



Xei(g),^^oi(g) 



^1 + il/V;[ya„(g)±K.a(q)] 



(16) 
t 

(17) 



«ci(g),Woi(g) = VfV^ + AIVc[Vaciq) ± K.a(9)(18) 

and Jireb(g) = Kob{q) = 1, fcb('7) = Wob(g) = vp for 6 > 1. 
Here we introduced a measure of the Coulomb interaction 
strength and defined the dimensionless parameter Vc — 
/ttvf£- 



C. Canonical transformation 

To treat the dot-lead interaction we use the canoni- 
cal transformation methods of Refs. 0S Following the 
approach of these papers but expressing the dot-lead in- 
teraction in terms of even and odd boson fields yields 

i^Coul = /^^^ E *9('^cl(<7) [WUq) + WRiq)] 



<l>ol{q)[WL{q)-WR{q)] ] 



(19) 



with Wa{q) = \og{l + l/ d^q'^). We observe that for a dot 
placed symmetrically between the two leads, the coupling 
to the odd mode vanishes exactly. 

The linear coupling between dot occupancy rid and 
lead density may now be removed by a canonical 
transformationi^ which shifts 



4>aiiq) 4>ai{q) - Za{q)nd/iqV2M 



with 



Zc{q) = ^^^^^^^MV, [WL{q) + WnM 



Zo{q) 



■Vc\{q) 
Voi{q) 



MVc[WL{q)-WBM 



(20) 

(21) 
(22) 



The canonical transformation acts on an operator O by 
O 



e'^Oe-'^ with 



S 



-na- 



/2ML 



EE^"('^)^H^- (23) 



iq 



-c.o q 

Under the canonical transformation the dot operator be- 
comes 

4 ^ 4e7ife5:,i^[nci(9)z„(g)+n„i(9)2o(9)]^ ^24) 



III. ELECTRONIC TUNNELING 

The crucial quantities in our considerations are expec- 
tation values of the form 



Fo.{t) = {c^pmo^xm) (25) 

with Yliaxp Ta^a\i3 + h-C- the renormalized hybridization 
Hamiltonian e^^ H^\y^e^^ without the da operator. We 
obtain, say, for a = L 

i^L(t) =i^o(t)e*^(*)e-*''^* (26) 

with Ffjit) the frce-fermion correlation and 



^ , , 27r v-^ 1 



-Sc(<z) 1 



-iVai{q)\q\t 



UF\q\t 



Bo{q)(l 



-iVai{q)\q\t 



(27) 



where we have defined 

Mq),B.iq)^\[K..iq)+^l^^ 



Koiiq) + 



[1 T Z,{q)Y 
Koliq) 



(28) 



To obtain an idea of the effects of screening, we ap- 
proximate the logarithmic functions Kei{q)-, Koi{q) and 
Vci{q), Voi{q) by constants Kd, i^oi and Vci, Woi- In this 
case the decay of electronic correlations is described by 



[-^sinhinTt)]- 



(29) 



with the Luttinger exponents 



1 / Kel 



(30) 

This formula plays a crucial role in our subsequent anal- 
ysis. We see that the odd-parity channel shift Zo (whose 
sign depends on which lead is nearer to the dot) enters 
with opposite sign in the two exponents. 

Equation (|5D|) implies a strong effect of the geometry 
on the electronic tunneling. In the symmetric case with 
equal dot-lead distances (11 = da = d/2, we see Zo = 0, 
because the odd boson mode does not contribute to the 
dot-lead coupling in that case. As the dot is moved off 
center, one of the exponents increases and the other de- 
creases; we shall see that this has important consequences 
for the tunneling. The dependence of the exponents on 
inter-lead distance and on the relative position of the lead 
is shown in Fig. [2] 

It is interesting to compare Eq. ([30| to the expressions 
obtained for the models studied in our previous workJ^ 
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(a) 




which is consistent with our result for the Luttinger expo- 
nent obtained in Ref. except that in this paper bound- 
ary rather than bulk exponents appeared. The crucial 
differences are the factor of 2 in the denominator and the 
factor of Z^/K which represented the 'orthogonality' ef- 
fect caused on one lead by adding a charge from the other. 
This orthogonality effect acts to suppress electron trans- 
port. By correlating the fluctuations in the two leads, 
the inter-lead Coulomb interaction acts to suppress this 
orthogonality effect; we shall see that the result is an 
enhanced tunneling. 

These arguments are supported by Fig. [2l The upper 
panel shows that the exponent increases with increas- 
ing inter-lead distance d, reflecting the weakening of the 
inter-lead interaction, the resulting enhancement of rel- 
ative fluctuations, and thus an increasing orthogonality 
effect. The lower panel shows that there is a relative 
dot-lead distance at which the exponent for tunneling 
from a given lead is minimized, reflecting the interplay 
between the increase in Zo and the decrease in Z^ due to 
the weakening of one of the couplings. For a very asym- 
metric situation we revert to the single-lead case studied 
previously. 



IV. CURRENT- VOLTAGE CHARACTERISTICS 



FIG. 2: (a) Luttinger exponent Yl of a symmetric junction 
[Yr — Yl) with nanotube leads as a function of the inter- 
lead distance d for different q values, (b) Luttinger exponents 
Yl, Yr as a function of the displacement 5 = di — d/2 of the 
quantum dot relative to the leads. [Here 5 = represents a 
symmetric junction. 



Reference^ considered a dot coupled to a single Luttinger 
liquid. In this case we have only one channel, and we 
found 



Y = — i M -1 
M \ 



K 



(31) 



The assumption of weak hybridization allows us to in- 
vestigate the nonequilibrium dynamics of the system us- 
ing a master-equation approach. The quantum dot with 
a spin-degenerate level is described by the diagonal den- 
sity matrix, which contains the occupation probabilities 
Po, Pi, and ■p2 of the empty state |0), the singly-charged 
state |1), and the doubly-charged state |2). Expand- 
ing the von Neumann equations to lowest order in the 
hybridization and making use of the Markov approxi- 
mation yields a set of master equations as discussed in 
Refs. 0@ with the tunneling rates related to the correla- 
tion function in Eq. (1251) . For example, the rate to move 
an electron from lead a to the empty dot is given by 
7?.Q_j.j = 2Re drFQ(r) e^*^'*'^. Evaluating the expres- 
sion we find that 



which can be obtained from our Eq. (j30p by dropping, 
say, all of the right channel interactions and couplings, 
which in practice means replacing Ze and Zo by Z , and 
Kei and K^i by K. The physics is therefore relatively 
similar between the two cases, except that because the 
odd-channel interactions are weaker than the even chan- 
nel ones, the exponent can be a little closer to the non- 
interacting value. 

However, if the quantum dot is coupled to two leads 
but we ignore the inter-lead electron-electron interac- 
tions, we may replace Zc by + Zfl, Zo by Zl — Z^ 
and /Cei, i^o\ by K to obtain 
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(32) 



\£d - ^J■o 

VrA 



(33) 



in the zero-temperature limit. Here 1/tq, cx 
27r|7^P/(uc/A) is the bare tunneling rate, and the Lut- 
tinger parameters have been approximated by constants. 

Solving the set of master equations in the steady state 
{Vi = 0) gives an expression for the steady-state current 
in terms of the transition rates. In the limit of very large 
U, the rates for the excitation of two electrons on the 
dot vanish and the current is described by the simple 
expression 



(/a) = 2e 
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FIG. 3: Comparison of current-voltage characteristics of 
a quantum dot symmetrically side-coupled to two parallel 
Luttinger liquid leads (a) for a model including inter-lead 
Coulomb interactions and (b) for the same model but with 
inter-lead Coulomb interactions set to zero. Voltages and 
thermal energies are given in units of the corresponding onsite 
energy. For the results shown in (a), the interaction parame- 
ters are obtained from Eq. (O using Eqs. ([T7} — (HH)) and 
with Vc 2i 0.9, Af = 4 (channels), da = d/2 = A, and with 
the logarithms evaluated at g = 0.1 (these parameters 
were shown in Ref. @] to be appropriate for a carbon nan- 
otube with A of the order of the tube diameter). The results 
shown in (b) are obtained for the same parameters except 
Vaa [cf. Eq. (|16|) ] set to zero. The local Coulomb interaction 
is assumed to be twice as large as the onsite energy in both 
panels. 



with TZn~^rn = X^a^n^m- Our numerical results are 
obtained with the full expression, which is presented in 
Ref. Also in our numerical results we used the param- 
eters identified in Ref. |^ as corresponding to a nanotubc 
with the Coulomb interaction cut off at a momentum 
scale of one tenth of the inverse tube diameter. 

The two panels of Fig. [3] compare the current- voltage 
characteristics of a quantum dot symmetrically side- 
coupled to two parallel Luttinger liquid leads with the 
inter-lead interactions included (upper panel) and ne- 
glected (lower panel). We see that including inter-lead 
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FIG. 4: (a) Current-voltage characteristics of a symmetric 
junction with nanotube leads for different inter-lead distances 
d. (b) Current-voltage characteristics for different displace- 
ments 5 = dL — d/2 of the quantum dot. [Here 5 = repre- 
sents a symmetric junction.] We assume the same parameter 
values as used in Fig. [S] 



Coulomb interactions in our model enhances the elec- 
tronic tunneling and changes the shape of the Coulomb- 
blockade steps qualitatively in comparison to the case of 
no inter-lead Coulomb interactions. The enhanced tun- 
neling, which manifests itself as regions of negative dif- 
ferential conductance, is due to the Coulombic dot-lead 
interaction, Eq. ([6]), as we discussed in Refs. and H. 
However, in the regime of weak electron-electron inter- 
actions, this enhancement can be overcompensated by 
orthogonality effects, which results in an effective sup- 
pression of the current. In general, the presence of inter- 
lead Coulomb interactions increases the total interaction 
strength and thus additionally enhances transport. The 
value of Ya is thus smaller than it would be if inter-lead 
Coulomb interactions were ignored. 

In the upper panel of Fig. |4]we demonstrate the effect 
of weakening the inter-lead interaction by moving the 
wires further apart. Increasing d weakens the strength 
of the interaction and thus suppresses the tunneling of 
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electrons. Increasing the inter-lead distance decouples 
the leads, causing the power-law 'divergences' explained 
above to vanish for sufficiently large inter-lead distances. 
Mathematically, the Luttinger exponents Ya then assume 
values larger than unity. 

Similarly, the lower panel of Fig. |4] demonstrates that 
breaking the spatial symmetry of the lead-dot-lead sys- 
tem can cause a suppression of the current near the 
Coulomb-blockade threshold, which would not occur in 
a quantum dot symmetrically side-coupled to the two 
leads. The dependence of the Luttinger exponent on 
the position of the dot relative to the leads shows that 
a displacement of ~ 10% already suffices to destroy the 
effect of enhanced electronic tunneling, since Yl then as- 
sumes values larger than unity. 

V. CONCLUSIONS 

We have studied transport through a quantum dot 
side-coupled to two parallel Luttinger liquid leads in the 
presence of a Coulombic dot-lead interaction. The new 
physics considered in this paper is the inter-lead inter- 
action. The geometry is chosen to allow a bosonisation 
solution to the inter-lead interaction as well as to the 
interaction between charge fluctuations on the dot and 
the dynamically generated image charge in the leads. A 
master-equation approach that treats the hybridization 



perturbatively has been applied to compute the tunnel- 
ing current. 

Our most important finding is that for symmetrical 
junctions the inter-lead electron-electron interactions en- 
hance the electronic transport and for reasonable in- 
teraction strengths may even change the overall sign 
of the Luttinger exponent Ya — I such that a suppres- 
sion (Ya > 1) of the current near the Coulomb-blockade 
threshold in the absence of inter-lead interactions turns 
into an enhancement (Ya < 1) in the I-V curves. The 
enhancement manifests itself as regions of negative dif- 
ferential conductance. As the symmetry of the junction 
is broken by moving the dot closer to one lead, the effect 
decreases. 

An important topic for future research is extending 
our results to the case of an end-coupled dot, where the 
inter-lead interaction breaks the mathematical transla- 
tion invariance which is employed in the standard theory 
of edge effects in Luttinger liquids. 
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